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Abstract 

ADM tetrad gravity is an Hamiltonian reformulation of General 
Relativity which gives new insight to the Dark Matter Problem. We 
impose constraints on the parameter space of ADM tetrad gravity with 
a Yukawa-like ansatz for the trace of the extrinsic curvature of the 3D 
hypersurfaces by htting the orbit of the S2 star around the Black Hole 
in the Galactic center and using the perihelia of some of the planets 
of the Solar System. We find very thight constraints on the strength 
of the coupling, 4.2 x 10”^ AU < 6 < 4.6 x 10“'^ AU, and an upper 
limit for the (inverse) scale length, fj, < 3.5 x 10~® AU“^. 


1 Introduction 


ADM tetrad gravity (ADM tg) [H [21 |3l HI [5l [6l [7] (first introduced in 
[8]) is an Hamiltonian reformulation of General Relativity based on a 
3-|-l splitting of spacetime into 3-dimensional spacelike hypersurfaces 
parametrized by the timejust like in the ADM formulation [9], but 
it considers tetrads as the dynamical variables instead of the metric. 

In [Q|7], authors studied Post Minkowskian (PM) and Post New¬ 
tonian (PN) expansion, in particular they calculated the equation of 
motion for a particle of mass m subjected to the gravitational force 
due to a potential <h, showing that there is a 0.5 PN (i.e. at the order 
0{c~^)) correction depending on the trace of the extrinsic curvature 
and on the velocity of the particle: 
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^These hypersurfaces are simultaneity surfaces. 
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where: 


( 2 ) 



with being the first order trace of the extrinsic curvature. 

Rearranging equation 0 , we have: 



(3) 


showing that the 0.5 PN correction can be interpreted also as an ad¬ 
ditional (position and velocity dependent) mass, thus giving an inter¬ 
esting and new insight on the Dark Matter problem, which, in this 
formalism, is an inertial effect due to the shape of the simultaneity 


hypersurfaces ^ * 


As has been proved in |^, in the ADM tg formalism, is a gauge 
variable, but there is no known natural gauge for this function: in this 
paper, we first discuss our ansatz for ^/C(i). 

In section]^ we fit equation Q to the data of the star S2 orbiting 
the black hole in the center of the Galaxy described in [min]; we will 
then find the part of the parameter space in which ADM tg gives a 
better fit than Newtonian gravity: this will impose limits on the free 
parameters of the theory (this is similar to what |14] did with f{R, (p) 
theories with a Sanders-like potential, m)- 

In order to improve our constraints, in section]^ we derive a for¬ 
mula for the precession angle given by the 0.5 PN correction and use 
it and perihelia precession angles of some of the planets of the Solar 
System to impose upper limits to our free parameters. 

Finally, in section we combine those constraints. 
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As we said in the introduction, in [3], it was shown that in ADM tg 
the trace of the extrinsic curvature is a gauge variable, but no natural 
gauge is known; we shall now suggest a possible ansatz, at least for 
the first order 

In [6l [7], authors also calculate the PN expansion of the metric, 
showing that the time-time component is given by (<I> is the Newtonian 

^As has been pointed out in 013, in this way, the Dark Matter becomes a metrology 
problem, since the trace and the shape of the hypersurfaces are linked to the choice of 
the convention on clock synchronizations. 

^We note further, that equation ([^ implies a violation of the Equivalence principle in 
the 3-space {not in the four dimensional spacetime). 
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potential): 


$ 2d 




( 4 ) 


It is known (see the review [TO] and references therein), that the 
PN expansion of /(i?) theories has the form: 


= 1 - 2 


GM 1 1 


r 1 + 5 


1 + (5 exp {—fJ-r) 


where the constant 5 is the measure of the Yukawa coupling. 

Confronting the previous equations one can make the ansatz that 
the trace of the extrinsic curvature has the form: 

= A{ct) - exp(—//r) 

and one can assume that the time dependence is linear A(ct) = ct 6, 
so our ansatz is: 

o- 1 

= cth - exp(—/rr). (5) 

In this way @ becomes: 

= 1 - 2 ^ - 2h ^Y(i)(ct,f). (6) 

There are two free parameters in ([^: the strength of the Yukawa 
coupling, 5, and the (inverse) length scale, /u; the aim of this work is 
to impose constraints on these parameters. 


3 Constraints from the S2 star 

In this section, we use publicly available data (see the complementary 
material of [T2|) of the orbit of the S2 star around the black hole in 
the center of the Galaxy to reduce the parameter space, in a similar 
fashion to m- 

3.1 Method 

We fit the orbit given by the equation 0 to the data reported in 
mm, using the ansatz discussed in the previous section (eqn. @)- 
We use the following algorithm, written with Wolfram Mathemat- 
ica 8, to carry out our fits: 

1. We start with an educated guess for the initial position of the 
star (the initial position and velocity of the star); 
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2. With those data we calculate the orbit of the star numerically 
with a 4*^ order Runge Kutta, obtaining the theoretical position 
of the star: xth-,- ■ ■; 

3. We minimize the (see HU): 

/ Xj - xth{i) \^ ^ ( Vi - ythii) ^ 

\ dxi ) V dyi J 

where Xi and yi are the observed positions and dxi and dyi are 
the relative errors; the index i runs over the obsevations. The 
parameters used in the fit are the initial conditions of the star 
and the strength of Yukawa coupling, 6, and it scale length, 

4. Errors are calculated using the Fisher Matrix method (see also 
[TT] and [T3]l. 

In our simulations we followed mm and in the Newtonian po¬ 
tential we fixed the mass of the central black hole (supposed fixed and 
point like) to = 4.3 x 10® M©; we also fixed the distance of the 
Sun from the center of the Galaxy to Rq = 8.3 kpc. 

The results of our fits are given in the next section. 

3.2 Results 

As a result of our fits, we obtain = 1.5477 for a Keplerian orbit 
and = 1.516890, with two additional degrees of freedom, for ADM 
tg: so the latter gives only a marginally better fit to the orbit. 

Orbits are shown in figure the best fit values for S and /j. are 
given in table together with their errors. 

In figure we plot the parameter region in which ADM tg gives 
a better fit to the orbit than Newton’s gravity: in green we plot the 
region where = 1.5477, the white dot is our best fit. 

We see that the parameters are not well constrained: we only have 
that (5 > 4.0 X 10“^ AU (see figure [l(b)] ) . 

In the section we shall combine these constraints with other 
coming from the precession or perihelia of some of the planets of the 
Solar System: this will leave available only a small region of parameter 
space. 



4 Constraints from the Solar System 

In HSl ES EH] the precession of perihelia is used to set constraints 
on the density of Dark Matter in the Solar System: the correction to 
the precession coming from Dark Matter must be at most of the same 
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Keplerian orbit 


1.5477 

ADM tg orbit 

Parameter 

Value 

Error 

5 (AU) 

0.00045 

0.0010 

p (AU)-i 

8.7 X 10”^ 

0.04 


1.516890 


Table 1: Fit results: the reduced th® Keplerian orbit and for ADM tg 

(we remind that the ADM tg has two more degrees of freedom). For ADM 
tg we give also the best £t values and the respective errors for the strength 
of the Yukawa coupling 6 and its legth scale p. 



(a) Light green: region of parameter 
space where ADM tg gives a better fit 
to the orbit then Newton gravity. 


(b) A zoom in the region: 3 x 10 < 

(5(AU) < 5.5 X 10-4 and 10"® < 

^(AU"^) < 10"^ . 


Figure 1: In both hgures: on the abscissa, the 5 on the ordinate /i. Green 
line = 1.5477 (our reduced Newtonian gravity). 


order of magnitude as the error on the measurements, thus imposing 
un upper value to its density; here we do the sameQ the correction to 
precession angle given by our 0.5 PN term (equation 0 , in the next 
section) will imopose upper limits on the couple 

We shall now derive a formula for the precession angle given by 
our 0.5 PN term. 

^After all we are aiming to describe Dark Matter as a manifestation of the non Eu- 
clidicity of the hypersurface. 
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(a) Keplerian orbit (b) ADM tg orbit 

Figure 2: Fit to the orbit of the S2 star: htted data and relative errors are 
in blue, the red line is the orbit. 

4.1 Precession angle 

As can be seen from equations Q, our perturbing force depends on 
the velocity x of the star, so it is not a central force. We follow m 
and use the Hamilton vector, since this mehod is easily generalized 
to our case (see [ISlESlISniEIlES] for the definition of the Hamilton 
vector and in particular, see the references [iziEaisi] for its use in 
the calculation of the perihelion precession angle). 

The Hamilton vector, u, is another conserved vector of the Kepler 
problem; it not independent from the more known Runge-Lenz R 
vector, in fact they are connected by the relation [17j0 

R = u X L 


where L is the angular momentum. 

Given a particle of mass m subjected to the gravitational force of 
a central object of mass M, the Hamilton vector and its magnitude 

^In this equation it is evident that the Hamilton vector is a conserved quantity, since 
it is linked to other constants of the motion. 
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are given respectively by HZ!: 


GM 


u = X — 


= r Cr 


IT"" 


r ip 


(7a) 


GM\ 


h ) 


, , GM 
ki = —-— e 
h 


(7b) 


where h is the magnitude of the angular momentum per unit mass 
and e is the eccentricity of the orbit and eV and are unit vectors {ip 
is the polar coordinate and r is the radial one). 

The precession rate for the Hamilton vector is given by [IZll25|: 


U X u 


w = 


\u\ 


where, (using the first order equation of motion): 




^ (^er + np ePj 


d? 


“ = ■ c ^ e- + r j /C(i) 


so, we find: 


w = — 


h r dj^ 


GMe c dt"^ 


( 8 ) 


where k is the unit vector orthogonal to the plane of the orbit. 
The precession is the integral of the rate ([8|) over a period T: 


A I I j ^ 

A9p = \w\dt = — 

Jo 
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. <iTi 


cGMe I 


(9) 


where r is the first order radial velocity (see [26]). The last equation 
can be reformulated as an integral over the radial coordinate r: 


Aep = -- ^ 


I 

Jr 


r+ ^2 

—TT dr 


c GMe J^_ dt'^ 

Now we can introduce the variable z [HllZllle]: 


( 10 ) 


5 


= 1 + e cos p = 1 + ez 


where S = h'^/GM is the semilatus rectum [Mini. With this variable, 
we have: 




c GM 




_s _ 

+ ez)^ \dt'^ 


dz 


( 11 ) 
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Planet 

Precession 
(xlO“^^ rad) 

Error 

(xl0“^^ rad) 

Mercury 

-1.5 

3.5 

Venus 

7.8 

4.8 

Earth 

0.9 

0.9 

Mars 

-0.19 

0.38 

Jupiter 

3400 

1632 

Saturn 

-45 

67 


Table 2: Precessions of planets per orbital period and their errors (see 

Eg). 


where: 





{dt^ + 2rdr + r^ ^^(i) = 




dr 


dt" + 2f — + — d 




dr 


2 + , ^2 + 2 


dd-2 


eS 


eS 


%) = 

d 


( 12 ) 


^JC 


( 1 ) 


Finally, from [26], we find r = h Vl — z'^. 


4.2 Results 

We use the data reported in table taken from |271128| : the error in 
Jupiter measurement is orders of magnitude bigger than the others 
and therefore it shall not be considered. 

Our results are given in figure]^ the black line is the upper limit 
set by Mercury, the green one by Venus, the blue one by the Earth; 
Mars and Saturn are almost overlapping (brown line); the allowed 
parameter space is colored in light blue. 


5 Combining both methods 

We can now combine the constriaints found in the previous sections. 

Referring to hgure [^ the black line is the combined upper limit 
imposed by Saturn and Mars, the light blue area is the same as in 
hgure the light green stripe is the area allowed by the orbit ht, 
hnally, the darker green area is the combined constraint. In formulae 
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Figure 3: Constraints on parameters imposed by planets’ precessions: on 
the abscissa 6 on the ordinate fi] black line is Mercury, green Venus, blue 
Earth; Mars and Saturn constraints are almost overlapping (brown line). 
The allowed parameter space is colored in light blue. 


we have: 



(13) 


6 Conclusion 

ADM tetrad gravity is an Hamiltonian reformulation of General Rel¬ 
ativity that gives new insight on Dark Matter problem: Dark Matter 
can be seen as an inertial effect due to the choice of the 3-1-1 splitting 
of the spacetime. 

In this work we suggested a possible ansatz for the first order 
of the trace of the extrinsic curvature: a Yukawa-like force with a 
linear time dependence. In this ansatz there are two free parameters: 
the strength of the Yukawa coupling 5 and its length scale. We then 
imposed constraints on the values of these parameters. 

In order to constraint the parameter space, we fitted the orbit of 
the S2 star around the black hole in the galactic center to the one 
given by ADM tg equations of motion and confronted the result with 
a pure Keplerian orbit (see also m)- 

We find that the ADM tg fit is only marginally better than that 
given by the simple Newtonian potential, in fact the reduced are 
1.5477 for the Keplerian orbit and 1.516890 for the ADM tg with two 
additional degrees of freedom. 
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Figure 4: Combination of constraints: ligth green stripe is the constraint 
from the S2 star, light blue area is the allowed parameter space from the 
Solar System, darker green is the combined constraint. The red dot is the 
best £t for the S2 orbit. 

We then used the perihelia of some of the planets of the Solar 
System to impose upper limits on 6 and p. 

Both the previous methods give unsatisfactory results, but their 
combination sets very tight limits on the strength of the Yukawa cou¬ 
pling, 4.2 X 10“^ AU < 5 < 4.6 x 10“^ AU, and an upper limit for 
the (inverse) scale length, p < 3.5 x 10“® AU“^. 
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